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After a brief review of the theoretical description of nuclei based on nonrelativistic 
many-body theory and realistic hamiltonians, these lectures focus on its application to 
the analysis of the electroweak response. Special emphasis is given to electron-nucleus 
scattering, whose experimental study has provided a wealth of information on nuclear 
structure and dynamics, exposing the limitations of the shell model. The extension of the 
formalism to the case of neutrino-nucleus interactions, whose quantitative understanding 
is required to reduce the systematic uncertainty of neutrino oscillation experiments, is 
also discussed. 



1. Introduction 

Over the past four decades, electron scattering has provided a wealth of information 
on nuclear structure and dynamics. Form factors and charge distributions have been 
extracted from elastic scattering data, while inelastic measurements have allowed 
for a systematic study of the dynamic response over a broad range of momentum 
and energy transfer.^ In addition, with the advent of the last generation of contin- 
uous beam accelerators, a number of exclusive processes have been analyzed with 
unprecedented precision (for recent reviews, see, e.g. Ref.El). 

In electron scattering experiments the nucleus il mostly seen as a target. Study- 
ing its interactions with the probe, whose properties are completely specified, one 
obtains information on the unknown features of its internal structure. 

The emerging picture clearly shows that the nuclear shell model fails to provide 
a fully quantitative description of the existing data. More realistic many body 
approaches, in which correlation effetcs are explicitely taken into account, appear 
to be needed to explain electron scattering observables. 

In neutrino oscillation experiments, on the other hand, nuclear interactions are 
exploited to detect the beam particles, whose kinematics is largely unknown. Using 
the nucleus as a detector obviously requires that its response to neutrino interac- 
tions be under control. Fulfillment of this prerequisite is in fact critical to keep the 
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systematic uncertainty associated with the reconstruction of the neutrino kinemat- 
ics to an acceptable levelP 

These lectures are aimed at providing an introduction to nuclear many-body 
theory and its application to the calculation of the electron- and neutrino-nucleus 
cross section. Section[2]contains a short overview of the underlying dynamical model 
and formalism. For pedagogical purposes, the nuclear response is first analyzed in 
the case of a scalar probe in Section [31 while the generalization to electron and 
charged current neutrino interactions is discussed in Sections |4] and EJ respectively. 
Finally, Section [5] is devoted to a summary of the main results. 

2. Non relativistic Nuclear Many-Body Theory 

In nuclear many-body theory (NMBT) the nucleus is viewed as a collection of A 
pointlike protons and neutrons, whose dynamics are described by the nonrelativistic 
hamiltonian 



where p,; and m denote the momentum of the i-th nucleon and the nucleon mass, 
respectively. The two body potential Vij is determined by fitting dcutcron properties 
and ~ 4000 precisely measured nucleon-nucleon (NN) scattering phase shitfs.Sl It 
turns out to be strongly spin-isospin dependent and non central, and reduces to the 
Yukawa one pion exchange potential at large separation distance. The inclusion of 
the three-nucleon interaction, satisfying ( Vijk ) <C {vij ), is required to account 
for the binding energy of the three-nucleon systems.^ 

The many body Schrodinger equation associated with the hamiltonian of Eq. ([T]) 
can be solved exactly, using stochastic methods, for nuclei with mass number A < 
12. The resulting energies of the ground and low-lying excited states are in excellent 
agreement with experimental data.^ 

It has to be emphasized that the dynamics of NMBT is fully determined by the 
properties of exactly solvable system, and does not suffer from the uncertainties 
involved in many-body calculations, which unavoidably make use of approxima- 
tions. Once the nuclear hamiltonian is determined, calculations of the properties of 
a variety of nuclear systems, ranging from deuteron to neutron stars, can be carried 
out without making use of any adjustable parameters. 

The main difficulty associated with the use of the hamiltonian of Eq.(IT|) in 
a many-body calculation lies in the strong repulsive core of the NN force, which 
cannot be handled within standard perturbation theory. 

In the shell model, this problem is circumvented replacing the interaction terms 
in Eq.([T|) with a mean field, according to 
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Within this scheme, the many-body Schrodinger equation reduces to a trivial single 
particle problem, and the ground state wave function can be written in the form 

%il,...,A)^A H , (3) 

aie{F} 

where the antisymetrization operator takes into account Pauli principle and the 
4>ai (i) are solutions of the eigenvalue equations 

(^ll + C/,^ 0a.W =ea,0a.W . (4) 

The product appearing in Eq.(l3]) includes the states belonging to the Fermi sea 
{F}, i.e. the A lowest energy states. 

In the Fermi gas (FG) model, in which neglecting interactions are neglected 
altogether, Eq.Q describes a degenerate gas of nucleons occupying all momentum 
eigenstates belonging to the eigenvalues |k| < kp. The Fermi momentum kp is 
related to the nucleon density p = A/V, V being the normalization volume, through 
kp = 2p/37r2. 

In spite of its simplicity, the shell model provides a remarkably good description 
of a number of nuclear properties. However, the results of electron- and hadron- 
induced nucleon knock-out experiments have provided overwhelming evidence of its 
inadequacy to account for the full complexity of nuclear dynamics. 

While the spectroscopic lines corresponding to knock-out from shell model orbits 
can be clearly identified in the measured energy spectra, the corresponding strengths 
turn out to be consistently and sizably lower than expected, regardless of the nuclear 
mass number. 

This discrepancy is mainly to be ascribed to the effect of dynamical correlations 
induced by the NN forces, whose effect in not taken into account in the shell model. 
Correlations give rise to virtual scattering processes, leading to the excitation of 
the participating nucleons to states of energy larger than the Fermi energy, thus 
depleting the shell model states within the Fermi sea. 

The first realistic treatment of the nuclear many-body problem was based on G- 
matrix perturbation theory, developed by Briickner, Bethe and Goldstone.^ Within 
this approach the nuclear hamiltonian is rewritten in the form (we neglect the three- 
body potential, for simplicity) 

H = Ho + Hi (5) 

with 

i ^ j>i i 

The single particle potential Wi is chosen in such a way as to either simplify the 
calculations or satisfy analitycity properties. The problem of the divergences arising 
from the short range repulsion of the bare interaction is circumvented replacing the 
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NN potential v with the G-matrix, obtained summing up ladder diagrams at all 
orders in the perturbative expansion through the integral equation 

G = v-v^ G, (7) 

where Q is a projection operator accounting for the effects of Pauli blocking and e 
is an energy denominator. 

A widely employed alternative approach exploits the possibility of embodying 
non perturbative effects in the basis functio ns. Th is is the foundation of Correlated 
Basis Function (CBF) perturbation theory, SMin! which provides a consistent and 
unified treatment of light nuclei and nuclear matter. 

The correlated states are obtained from the eigenstates of the shell model hamil- 
tonian through the transformation 

_ F\n) 
'""^ " (n|FtF|n)i/2 ' 

where -F is a correlation operator, whose structure reflects the properties of the NN 
potential. It is written in the form 

F = SY[f,,, (9) 

j>i 

where S is the symmetryzation operator and 

n 

The minimal set of operators Ofj includes the four central components associated 
with the different spin-isospin channels {n = 1,4) and the isoscalar and isovector 
tensor components (n — 5,6): 

Of^. = [1,(^,.^,),5,,]®[1,(t,.t,)] , (11) 

with 

3 

= ■ ~ ('^^ ■ '^j) ■ (12) 

In some calculation the spin orbit (L • S) and isovector spin orbit ((L • S){Ti ■ Tj)) 
components, needed to reproduce the NN scattering phase shifts in S and P wave, 
have been also included. 

The radial functions /" (r^ ) are determined through functional minimization of 
the expectation value 

(o™ 

evaluated using the cluster expansion formalism.!^ 

The correlated wave functions \n) are required to exhibit two important 
properties!^ the cluster factorization property, dictated by the finite range of the 
interaction, and the core property, due to the short range repulsion. The former 
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property implies that /"^^(r) 1 and /"^^(r) ^ as r cx), while the core 
property requires that all correlation functions become vanishingly small at short 
interparticle distance. 

The basis of correlated states, while being complete, is not orthogonal. However, 
it can be orthogonalized using standard techniques of many-body theory.^l^ As the 
operator HF is well behaved, as is the G-niatrix, the perturbative expansions of 
nuclear observable in the correlated basis are rapidly convergent. 

3. The nuclear response to a scalar probe 

Within NMBT, the nuclear response to a scalar probe delivering momentum q and 
energy w can be written in te rms of the the imaginary part of the polarization 
propagator n(q, w) according tJHE^ 

5(q, u;)^-lm n(q, c.) = ^ Im (0|pt ^ pq|0> , (14) 

TT TT ^H — Ff)— uo — irj 

where = 0^, pq = Sk'^k+q'^k is the operator describing the fluctuation of the 
target density induced by the interaction with the probe, aj^ and at are nucleon 
creation and annihilation operators, and |0) is the target ground state, satisfying 
the Schrodinger equation H\Q) = Fq\Q). 

In this Section, we will discuss the relation between ^(q, w) and the nucleon 
Green functio n, lead ing to the popular expression of the response in terms of spec- 
tral functions .fl2ll3! This discussion is mainly aimed at showing that the spectral 
function formalism, while being often advocated using heuristic arguments, can be 
derived in a rigorous and fully consistent fashion. For the sake of simplicity, we will 
consider uniform nuclear matter with equal numbers of protons and neutrons. 

Equation (fT4|) clearly shows that the interaction with the probe leads to a tran- 
sition of the struck nucleon from a hole state of momentum k, with |k| < kp^ to a 
particle state of momentum k 4- q, with |k 4- q| > kp. The calculation of ^(q, w) 
amounts to describing the propagation of the resulting particle-hole pair through 
the nuclear medium. 

The fundamental quantity involved in the theoretical treatment of many-body 
systems is the Green function, i.e. the quantum mechanical ar npl itude associated 
with the propagation of a particle from x = (t, x) to x' = (t', x') .Li^In uniform mat- 
ter, due to translation invariance, the Green function only depends on the difference 
X — x' , and after Fourier transformation to the conjugate variable k = ( k , _E ) 
can be written in the form 

G{k,E) = {0\al ^ ^ \ . ak|0) - (0|ak „ r, \ ■ «klO) 

H — Fo — F — 17] H — Fq + F — If] 

= Gh{k,F) + Gp{k,E) , (15) 

where Gh and Gp correspond to propagation of nucleons in hole and particle states, 
respectively. 
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The connection between Green function and spectral functions is established 
through the Lehman representatiorffH 



E' -E-irj E -E' -iT] 



(16) 



G(k, E)^ J dE' 
implying 

PhiKE) - J2 \{niN-i)i-k)\a^\ON)fS{E ~ E^^ + ^o) = ^ Im G„(k, E) (17) 

n 

Pp(k,£;)-^|(n(jv+i)(k)|a[|0Ar)|25(i? + i?(+)-i?o) = ^ImGp(p,£;) , (18) 

n 

where l(n(jv±i)(±k)} denotes an eigenstate of the (A± l)-nucleon system, carrying 
momentum ±k and energy eI^\ 

Within the FG model the matrix elements of the creation and annihilation 
operators reduce to step functions, and the Green function takes a very simple 
form. For example, for hole states we fincfl 

G^G,.(k,i?)-^^^o_.^ , (19) 

with e". = |k^|/2TO, implying 

PFG,h{^, E) = e{kF - \kmE + el) . (20) 

Strong interactions modify the energy of a nucleon carrying momentum k according 
to el — > el + S(k, E), where S(k, E) is the complex nucleon self-energy, describing 
the effect of nuclear dynamics. As a consequence, the Green function for hole states 
becomes 

1 

E + el+'Y.{k,E) ■ 



Gh{KE)= ^ , „ , . (21) 



A very convenient decomposition of G;i(k, E) can be obtained inserting a com- 
plete set of (A— l)-nucleon states (see Egs . (flSl) - () 1 8 p J and isolating the contributions 
of one-hole bound states, whose weight is given bjU4l 

Zk = |(-k|ak|0)|2 = d{kF - |k|)$fe . (22) 

Note that in the FG model these are the only nonvanishing terms, and <f>/j = 1, 
while in the presence of interactions <i>/c < 1. The resulting contribution to the 
Green function exhibits a pole at — efc, the quasiparticle energy Cfe being defined by 
the equation 

efe = e°+ReI](k,efe) . (23) 



'^Note that, according to our definitions, the hole spectral function is defined for E > — /i, /i being 
the Fermi energy. 
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The full Green function can be rewritten 

Gh (k, E) = J\ + Gl (k, E) , (24) 

+ Efe + zZfe Im S(k, Cfc) 

where is a smooth contribution, asociated with [A — l)-nucleon states having 
at least one nucleon excited to the continuum (two hole-one particle, three hole-two 
particles . . . ) due to virtual scattering processes induced by nucleon-nucleon (NN) 
interactions. The corresponding spectral function is 

' ^ 1^ + + Re Et!:^":;!^ e,.. + ^' ■ (^^) 

The first term in the right hand side of the above equation yields the spectrum 
of a system of independent quasiparticles, carrying momenta |k| < /cf, moving in 
a complex mean field whose real and imaginary parts determine the quasiparticle 
effective mass and lifetime, respectively. The presence of the second term is a con- 
sequence of nucleon-nucleon correlations, not taken into account in the mean field 
picture. Being the only one surviving at |k| > kp^ in the FG model this correlation 
term vanishes. 
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Fig. 1. Energy dependence of the hole spectral function of nuclear matter.HSl The solid, dashed 
and dot-dash lines correspond to |k| = 1, 0.5 and 1.5 fm~^, respectively. The FG spectral function 
at |k| = 1 and 0.5 fm~^ is shown for comparison. The quasiparticle strengths of Ea.l l22ll , are also 
reported. 



Figure [T] illustrates the energy dependence of the hole spectral function of nu- 
clear matter, calculated in Ref.'^ using CBF perturbation theory and a realistic 
nuclear hamiltonian. Comparison with the FG model clearly shows that the ef- 
fects of nuclear dynamics and NN correlations are large, resulting in a shift of 
the quasiparticle peaks, whose finite width becomes large for deeply-bound states 
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with |k| <^ kp. In addition, NN correlations are responsible for the appearance of 
strength at |k| > kp- The energy integral 

n(fc) = JdEPh{k,E) (26) 

yields the occupation probability of the state of momentum k. The results of Fig. 
[T] clearly show that in presence of correlations n(|k| > kp) 0. 

In general, the calculation of the response req uires th e knowledge of Ph and Pp, 
as well as of the particle-hole effective interactionP^HSl xhe spectral functions are 
mostly affected by short range NN correlations (see Fig. [T|) , while the inclusion of 
the effective in tera ction, e.g. within the framework of the Random Phase Approx- 
imation (RPA),liSl is needed to account for collective excitations induced by long 
range correlations, involving more than two nucleons. 

At large momentum transfer, as the space resolution of the probe becomes small 
compared to the average NN separation distance, S'(q, uj) is no longer significantly 
affected by long range correlations. In this kinematical regime the zero-th order 
approximation in the effective interaction, is expected to be applicable. Whithin 
this scenario, the response reduces to the incoherent sum of contributions coming 
from scattering processes involving a single nucleon, and can be written in the 
simple form 

S{ci,uj)^ J d^kdE Ph{k,E)Pp{k + q,uj-E) . (27) 

The widely employed impulse approximation (lA) can be readily obtained from the 
above definition replacing Pp with the FG result, which amounts to disregarding 
final state interactions (FSI) betwen the struck nucleon and the spectator particles. 
The resulting expression reads 

5/A(q,w) = J d^'kdE Ph{k,E)e{\k + q\-kp)Siu;-E~e^^^^^^) . (28) 

Figure [2l showing the lu dependence of the nuclear matter structure function at 
|q| = 5 fm^^, illustrates the role of correlations in the target ground state. The solid 
and dashed lines have been obtained from Eg. ([28]) using the spectral function of 
Ref.El and that resulting from the FG model (shifted in such a way as to account for 
nuclear matter binding energy), respectively. It clearly appears that the inclusion 
of correlations produces a significant shift of the strength towards larger values of 
energy transfer. 

At moderate momentum transfer, both the full response and the particle and 
hole spectral functions can be obtained using non relativistic many-body theory. 
The results of Ref."*^^ suggest that the zero-th order approximations of Eqs. ([27|) and 
([28|) are fairly accurate at |q| ^ 500 MeV. However, in this kinematical regime 
the motion of the struck nucleon in the final state can no longer be described using 
the non relativistic formalism. While at lA level this problem can be easily circum- 
vented replacing the non relativistic kinetic energy with its relativistic counterpart. 
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Fig. 2. Nuclear matter S'jA(q, toi) (see Ea. l|28| l'). as a function of ui at |q| = 5 fm ^. The solid and 
dashed lines correspond to the spectral function of Rem and to the FG model (shifted in such 
a way as to account for nuclear matter binding energy), respsctively. 



including the effects of FSI in the response of Eq. ([77]) involves further approxima- 
tions, needed to obtain the particle spectral function at large |q|. 

A systematic scheme to include corrections to Eq. (l^5|) and take into account 
FSI, originally proposed in Ref.lIS!^ jg discussed in Ref."'^'''. The main effects of FSI 
on the response are i) a shift in energy, due to the mean field of the spectator 
nucleons and ii) a redistributions of the strength, due to the coupling of the one 
particle-one hole final state to n particle-n hole final states. 

In the simplest implementation of the approach of Refsff^^^, the response is 
obtained from the lA result according to 



Sici,uj) = J duj' SiAici,uj')fq{u 



J) 



(29) 



the folding function /q being related to the particle spectral function through 

Pp(k + q,c^-i?)-0(fcf.-|k + q|) /|k+q|(c^-i?-efk+q|) (30) 



ith 



|k+q| 



= ^kT^ 



m?. In the absence of FSI, /q shrinks to a (5-function 



and the lA result of Eq. ([28|) is recovered. 

Obvioulsy, at large q the calculation of Pp (k -I- q, w — E) cannot be carried out 
using a nuclear potential model. Hovever, it can be obtained form the measured 
NN scattering amplitude within the eikonal approximation. The resulting folding 
function is the Fourier transform of the Green function describing the propagation 
of the struck particle, travelling in the direction of the z-axis with constant velocity 
v. 



dt 
2^ 



^iwt i/g dt'V[q|(t)t') 



(31) 
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where k + q « q and 

l^q|W = (0|^Er|q|(r«. +^)|0> • (32) 

In the above equation, r|q| is the Fourier transform of the NN scattering ampli- 
tude at incident mometum |q| and momentum transfer |t|, A|q|(fc), parametrized 
according to 

A\^\{p) ^ ^-^<J{^ + a)e-^''' . (33) 

In principle, the total cross section a, the slope (3 and the ratio between the real 
and the imaginary part, a, can be extracted from NN scattering data. However, 
the modifications of the scattering amplitude due to the presence of the nuclear 
medium are known to be sizable, and must be taken into account. The cal culation 
of these corrections within the framework of NMBT is discussed in Ref.[lHI. 

In Eq. ([32|) . the expectation value is evaluated in the correlated ground state. 
It turns out that NN correlation, whose effect on Ph is illustrated in Fig. [U also 
affect the particle spectral function and, as a consequence, the folding function of 
Eq. pO|l . Neglecting all correlations 

V\^\iz)^m^ = ^vpa{z + a) , (34) 
and the quasiparticle approximation 

Pp{q,iu-E)^-- ^1-2 -2 (35) 



E-e°-Rc 

|q| |q| 



Im V° 
|q| 



is recovered. 

Correlations induce strong density fluctuations, preventing two nucleon from 
coming close to one another. The joint probability of finding two particles at posi- 
tions ri and r2 can be written 

p(ri,r2) =p^5(|ri -r2|) (36) 

where the g(r) is the radial distribution function, shown in FigO 

The effect of correlation on FSI can be easily understood: as the probability 
of finding a spectator within the range of the repulsive core of the NN force (< 1 
fm) is small, the probability that the struck particle rescatter against one of the 
spectators within a length ^ l/|q| (the space resolution of the probe) is also very 
small at large |qj. H ence, inclusion of correlations leads to a significant suppression 
of FSI effects.ll^ 

Fig. [4] shows the uj dependence of the nuclear matter response of Eq. ((29|) at 
|q| = 5 fm"-'^. The solid and dashed lines have been obtained using the spectral 
function of Rem, with and without inclusion of FSI according to the formalism 
of Rei^, respectively. For reference, the results of the FG model are also shown 
by the dot-dash line. The two effects of FSI, energy shift and redistribution of the 
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Fig. 3. Radial distribution function of nuclear matter at equilibrium density, obtained from CBF 
perturbation theory using a realistic hamiltonian. 
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Fig. 4. Nuclear matter ^{q, a;) , defined in Eq. I I29I I , as a function of lij at | q| = 5 fm~ ^ . The solid and 
dashed lines have been obtained from the spectral function of RefU^, with and without inclusion 
of FSI, respectively. The dot-dash line corresponds to the FG model (shifted in such a way as to 
account for nuclear matter binding energy), respectively. 

strength from the region of the peak to the tails, clearly show up in the comparison 
betweem soild and dashed lines. 



4. The electron-nucleus cross section 

The differential cross section of the process 

e + A-^ e + X , 



(37) 
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in which an electron of initial four-momentum = {Ee,^e) scatters off a nuclear 
target to a state of four-momentum k'^ = (£'p/,kp/), the target final state being 
undetected, can be written in Born approximation aA^I 

where a — 1/137 is the fine structure constant, dVle' is the differential solid angle 
in the direction specified by kg/, = —q^ and q = k^ — ke> = (i-iJ,q) is the four 
momentum transfer. 

The tensor L^t^, that can be written neglecting the lepton mass as 

L^, = 2 [fc'ffc.n + k'^k^, - g^^'ikek,,)] , (39) 

where g'^" = (1,-1,— 1,-1) and {keke') = E^Ee' — kg • kg/, is fully specified by the 
measured electron kinematical variables. All the information on target structure is 
contained in the tensor W^'^ , whose definition involves the initial and final nuclear 
states |0) and carrying four-momenta po a-nd px, as well as the nuclear current 
operator J'^: 

W^^ =Y,{Q\J^\X){X\r\Q)5'^^Hpa + q-px) , (40) 

X 

where the sum includes all hadronic final states. Note that the tensor of Eq. (l40|) is 
the generalization of the nuclear response, discussed in Section [3l to the case of a 
probe interacting with the target through a vector current. This can be easily seen 
inserting 

EHH = i, (41) 

n 

|n) being an eigenstate of the nuclear hamiltonian sasisfying H\n) = En\n), in the 
definition of Eq. (fT4|) . The result is 



S{ci,u) = Y.'^Q\p\\n){n\p^\0)5{u + Eo - E^) , (42) 

n 

to be compared to Ea. (|^ . 

The most general expression of the target tensor of Eq. flQ]) . fulfilling the re- 
quirements of Lorentz covariance, conservation of parity and gauge invariance, can 
be written in terms of two structure functions Wi and W2 as 

= H- itf ) + i|(ps - '^,') [a - ^,-) , («) 

where M is the target mass and the structure functions depend on the two scalars 
and (poq)- In the target rest frame (po?) = and Wi and W2 become functions 

of the measured momentum and energy transfer |q| and lo. 
Substitution of Eq. (gS]) into Eq. leads to 

d'^cr ( da 



d^p'dEpi V dQp' 



M 



2 



T4^2(|q|,w) + 2P^i(|ql,a;)tan- ^ 



(44) 
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where 9 and {da/dfle')M — c? cos^(0/2)/4i?g sin"'(6'/2) denote the electron scatter- 
ing angle and the Mott cross section, respectively. 

At moderate momentum transfer, typically (|q| < 0.5 GeV), the response 
tensor W^'^ of Eq. ((40|) can be obtained from NMBT, using nonrelativistic wave 
functions to describe the initial and final states and expanding the current operator 
in powers of |q|/m. 

Exact calculations of the electron-nucleus cross section can be carried out for 
light nuclei, with ^ < 4, either solving the Schrodinger equati on for bound and con- 
tinuum stateJ2nior using integral transform techniques I ^^ I ^^ ^The latter approach is 
discussed in detail in Prof. Leidemann's lectures .^221 Accurate calculations are also 
possible for uniform nuclear matter, as translation invariance considerably simpli- 
fies the problem. On the othe hand, the available results for medium-heavy 
targets have been mostly obtained using the mean field approach, supplemented 
by the inclusion of model residual interactions to take into account long range 
correlations.[2S! 

At higher values of |q|, corresponding to beam energies larger than '^l GeV, 
describing the final states \X) in terms of nonrelativistic nucleons is no longer 
possible. Due to the prohibitive difficulties involved in a fully consistent treatment 
of the relativistic nuclear many-body problem, calculations of W^^ in this regime 
require a set of simplifying assumptions, allowing one to take into account the 
relativistic motion of final state particles carrying momenta ~ q, as well as inelastic 
processes leading to the production of hadrons other than protons and neutrons. 

4.1. The impulse approximation 

Within the lA picture, schematically represented in Fig. [5l the nuclear current 
appearing in Eq. (j40|) is written as a sum of one-body currents 

J'^T.^'^ (45) 

i 

while \X) reduces to the direct product of the hadronic state produced at the elec- 
tromagnetic vertex, carrying momentum p^;, and the state describing the residual 
system, carrying momentum p-^ = q — Pa; (in order to simplify the notation, spin 
indices will be omitted) 

|X)^|2;,p,)®|7^,p7^) . (46) 
Using Eq. ((46l) we can replace 

Y,\X){X\^Y1 I rfVk,P.)(p.,a:|^dV|7^,PK)(P7^,7^| • (47) 
X X •' n 

Substitution of Eqs. (|45p - (|47|) into Eq. ((40|l and insertion of a complete set of free 
nucleon states leads to the factorization of the nuclear current matrix element. As 
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Fig. 5. Schematic representation of the lA regime, in which the nuclear cross section is replaced 
by the incoherent sum of cross sections describing scattering off individual nuclcons, the recoiling 
(A — l)-nucleon system acting as a spectator. 

a result, the incoherent contribution to Eq. (|40p can be rewritten in the form 

W^''ici,iu) = J d'k dE (^-^^ [ZP,{k,E)w^''{q)+NP^ik,E)wi^,''{q)] , (48) 

where Z and N — A — Z are the number of target protons and neutrons, while 
Pp and Pn denote the proton and neutron hole spectral functions, respectively. In 
Eq. (gHl), Ek = ^|k2| +m2 and 

= ^(k, N| |x, k + q) (k + q, x\f^\N, k)<5(5 + E^ ^ E^) , (49) 

X 

witlffl 

uj ^ Ea; - E]^ = Eo + uj - Etz - Eu = uj - E + m - Ey, . (50) 

The above equations show that within the lA scheme, the definition of the electron- 
nucleus cross section involves two elements: i) the tensor w'^ , defined by Eq. (j49p . 
describing the electromagnetic interactions of a bound nucleon carrying momentum 
k and ii) the spectral function, discussed in the previous Sections. 

4.2. Electron scattering off a bound nucleon 

While in electron-nucleon scattering in free space the struck particle is given the 
entire four momentum transfer q = (w, q), in a scattering process involving a bound 
nucleon a fraction Suj of the energy loss goes into the spectator system. This mech- 
anism emerges in a most natural fashion from the lA formalism. 

Assuming that the current operators are not modified by the nuclear environ- 
ment, the quantity defined by Eq. (|49|) can be identified with the tensor describing 
electron scattering off a free nucleon at four momentum transfer q = (q, lu). Hence, 
Eq. (j49|) shows that within lA binding is taken into account through the replace- 
ment 

q= {uj,q) -^q= (w, q) . (51) 
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The interpretation oi Slo — lu ^ lu as the amount of energy going into the recoiling 
spectator system becomes particularly transparent in the limit |k|/r7i <^ 1, in which 
Eq. dnni) yields 6uj = E. 

The tensor w'^ of Eq. ([49]) can be obtained from the general expression (compare 
to Eq. 



< - < + ) + S - f^r) {>■■ - fr) . (52) 

where k = (i?k, k) and the two structure functions wi and W2 can be extracte d from 
the measured electron-proton and clectron-deuteron scattering cross sectionJ^^^^. 

For example, in the case of quasielastic scattering wi and W2 are simply related 
to the electric and magnetic nucleon form factors, Ge^ ^^id Gm^, through 

—2 / ^ \ 



A similar expression can be used to describe resonance production. The explicit 
formula differ in the analytical form of the relevant form factors and for the replace- 
ment of the energy conserving (5-function with a Breit-Wigner factor, accounting for 
the finite width of the rcsonance.^^SI Finally, Eq. ([52]) can be applied in the region of 
deep inelastic scattering, where the structure functions — mwf and F2 — ujui^ 
depend on = |q^| — and the bjorken scaling variable x = j^muj. 

As a final remark, it has to be pointed out that the replacement of lo with a), 
while being reasonable on physics grounds, and in fact quite natural in the context 
of the lA analysis, poses a considerable conceptual problem, in that it leads to a 
violation of current conservation, that requires 

g^'^TV = • (55) 

However, violation of gauge invariance in the lA scheme turns out to be only 
marginally relevant to inclusive electron scattering at large momentum transfer. 
The results of numerical studies suggest that the main effect of nuclear binding can 
indeed be accounted for with the replacement w — ^ tD.tH 



4.3. Comparison to electron scattering data 

The formalisr n outli ned above, supplemented by the treatment of FSI effects pro- 
posed in Refs.l^^i^, has been widely and successfully applied to the analysis of 
electron-nucleus scattering data.^^ 

In Ref. it has been employed to calculate the inclusive electron scattering 
cross sections off oxygen at beam energies ranging between 700 and 1200 MeV and 
electron scattering angle 32°. In this kinematical region, relevant to many neutrino 
experiments, single nucleon knock out is the dominant reaction mechanism and 
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both quasi-elastic and inelastic processes, leading to the appearance of nucleon 
resonaces, must be taken into account. 

Comparison between theoretical results and the experimental data of Ref. ^201 
shows that, while the data in the region of the quasi-elastic peak are accounted for 
with an accuracy better than ~ 10 %, theory fails to explain the measured cross 
sections at larger electron energy loss, where A prod ucti on dominates. 

As an example, Fig. [S] shows the results of Ref.l^H at beam energy 700 and 
1200 MeV. For reference, the results of the Fermi gas (FG) model corresponding to 
Fermi momentum pf = 225 MeV and average removal energy e — 25 MeV are also 
shown. Theoretical calculations have been carried out using the spectral functi on 
of Ref.E2l, the Hohler-Brash parameterization of the nucleon form factors I '^'^ l '^^ l in 
the quasi-elastic channel and the Bodek and Ritchie parametrization of the proton 
and neutron structure functions in the inelastic channels 




Fig. 6. Cross section of the process ^®0(e,e') at scattering angle 32° and beam energy 700 MeV 
(left panel) and 1200 MeV (right panel), as a function of the electron energy loss ui. Solid lines: full 
calculation, including FSI. Dot-dash lines: lA calculation. Dashed lines: FG model withpi? = 225 
MeV and e = 25 MeV. The data are taken from Ref^ 



The authors of Ref. argued that the disagreement between theory and data in 
the A production region is likely to be imputable to deficiencies in the description 
of the neutron structu re func tions at low (5^0 This conclusion is supported by the 
analysis carried out in.EHESl 

The left panel of Fig. [7] shows that the neutron structure function extracted 
from Jefferson Lab data at ~ 0.5 GeV^'^, following the procedure of Bodek 
and Ritchie, is significatly larger than the one resulting from the analysis of Ref.^^, 
based on SLAG data spanning the kinematical domain 1 < < 20 GeV^. The 
theoretical cross section obtained using the neutron structure function of 

Ref.I2Hl, 



''In the kinematics of Fig.[6]the A production peak corresponds to 0.2 GeV'^. 
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displa yed in the right panel, turns out to be in close agreement with the SLAC data 
of Ref!^, corresponding to ^ 0.4, in the A region. 




¥ [GeV] u [GeV] 



Fig. 7. Left panel: nucleon structure functions W2 (N = n,p) at Ee = 2.445 GeV and de = 20°, 
corresponding to = 0.54 GeV^ at the A production peak, plotted as a function of the invariant 
mass of the hadronic final state. The shaded area represents the ui^ resulting from the analysis 
of Ref.ESlj while the dashed and dot-dash lines correspond to ^^'^ '"'2 °^ 

Ref.|2Zl, respectively. 

Right panel: electron scattering cross section oflF carbon at Ee = 1.3 GeV and 9e = 37.5°, as a 
function of the electron energy loss lo. The solid line correspo nds to theoretical calculations carried 
out using the proton and neutron stru cture functions of Ref . whil e th e shaded region has been 
obtained using «)" and liij of Ref. ESI. The data are taken from Ref. 



5. Charged current neutrino- nucleus interactions 

In Born approximation, the cross section of the weak charged current process 

i^i + A^r+X, (56) 
can be written in the form (compare to Eq. (|38p ) 

da G2 Ik'l 



dQfdEf 327r2 Ikl 



W.W^'" , (57) 



where G = GpcosOc, Gp and 9c being Fermi's coupling constant and Cabibbo's 
angle, is the energy of the final state lepton and k and k' are the neutrino and 
charged lepton momenta, respectively. Compared to the corresponding quantities 
appearing in Eq. p8|l , the tensors L^^ and W^^'^ include additional terms resulting 
from the presence of axial- vector components in the leptonic and hadronic currents 
(see, e.g., Ref. ^2%. 

Within the lA scheme, the cross section of Eq. ([57|) can be cast in a form similar 
to that obtained for the case of electron-nucleus scattering. Hence, its calculation 
requires the nuclear spectral function and the tensor describing the weak charged 
current interaction of a free nucleon, . In the case of quasi-elastic scattering. 
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neglecting the contribution associated with the pseudoscalar form factor Fp, the 
latter can be written in terms of the nucleon Dirac and Pauli form factors -Fi and 
F2, related to the measured electric and magnetic form factors Ge and Gm through 
1 / g2 \ I 
■^1 = ^ TTa — 2 - -. — ^Gm , ^2 = yTa — 2 (^*^ ^ '^■e) , (58) 

and the axial form factor Fa- 

It has to be pointed out that the formalism described in Section [3l while includ- 
ing dynamical correlations in the final state, does not take into account statistical 
correlations, leading to Pauli blocking of the phase space available to the knocked- 
out nucleon. 

A rather crude prescription to estimate the effect of Pauli blocking amounts to 
modifying the spectral function through the replacement 

P{p,E)^P{p,E)9{\p + q\-pp) (59) 

where pp is the average nuclear Fermi momentum, defined as 

Pf = J d\ pa{v)pp{y), (60) 

with ^^(r) = (37r^/9A(r)/2)^/'^, pa{y) being the nuclear density distribution. For 
oxygen, Eq. ([60)) yields pp — 209 MeV. Note that, unlike the spectral function, 
the quantity defined in Eq. ((59)) does not describe intrinsic properties of the target 
only, as it depends explicitely on momentum transfer. 

The effect of Pauli blocking is hardly visible in the energy loss spectra shown in 
Fig. m as the kinematical setup corresponds to > 0.2 GeV^ at the quasi-elastic 
peak. However, it becomes appreciable at lower . 




Fig. 8. Quasi-elastic differential cross section da/dEfj,, as a function of the charged lepton energy 
Eij,, for neutrino energy E = 0.8 (left panel) and 2.0 GeV (right panel). The solid and dot-dash 
lines show the results of the lA calculation with and without Pauli blocking (implemented as in 
Eq. l|59|ll. respsctively. The dashed lines have been obtained using the FG model. 
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As an example, Fig. [8] shows the I'^-nucleus cross sections, as a function of the 
scattered muon energy, cal cul ated within the FG model (dashed line) and using the 
the spectral functio of RefP, with and without Pauli blocking (solid and dot-dash 
lines, respectively). It clearly appears that the FG model yields a larger peak at 
high-energy. This feature should show up in the cross section at forward angles, 
and may have a direct effect on neutrino oscillation measurements. 

6. Summary 

The approach based on NMBT provides a unified parameter-free description of the 
electroweak nuclear response in a variety of kinematical regions. 

Thanks to the availability of reliable spectral functions, accurate calculations 
of the cross sections in the lA regime can be carried including the effects of short 
range NN correlations. Correlation effects can also be consistenly included in the 
treatment of FSI between the struck nucleon and the spectator particles. 

Comparison to electron- nucleus scattering data shows that, the region of quasi- 
elastic scattering and A production can be reproduced with remarkable accuracy. 
The role of meson exchange currents, which are known to provide a significant 
amount of strength in the dip region between the quasi-elastic and the A peak, still 
needs to be carefully investigated. The generalization of the formalism discussed in 
these lectures to deep inelastic scattering is straightforward.^ 

As a final remark, it has to be pointed out that the possibility of using the 
approach based on NMBT in the analysis of neutrino experiments largely depends 
on the ability to implement its elements in Monte Carlo simulations. 

Assuming, for the sake of simplicity, that the elementary weak interaction vertex 
in the nuclear medium be the same as in free space, a realistic simulation of neutrino- 
nucleus scattering requires the energy and momentum probability distribution of 
the nucleons, needed to specify the initial state, as well as their distribution in space 
and the medium modified hadronic cross section, needed for the description of FSI. 

Studies based on NMBT and stochastic methods to solve the many-body Schro- 
dinger equation appear to be capable of providing access to all the above quantities 
for a variety of nuclear targets. 
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